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^ Abstract 

^ We study the magnetic ordering of Vanadium spinels by Quantum Monte Carlo simulations of a 

^ three-band Hubbard model. Vanadium spinels, AV2O4, exhibit a unique "up-up-down-down" spin 



ordering at low temperatures. While this magnetic ordering was originally measured in 1973, its 



tin 

origin has remained unclear for many years due to the lack of unbiased approaches for solving the 



relevant model. A three-band Hubbard model on the spinel lattice (corner sharing tetrahedra) is a 
minimal Hamiltonian for describing the t2g electrons of the V^"^ ions. One of the main difficulties is 
j_J that this family of compounds belongs to the elusive intermediate-coupling regime (U ^t) for which 

^ there is no small parameter that can justify a perturbative expansion. We present a controlled 

^ Quantum Monte-Carlo approach to the three-band Hubbard model relevant for this materials that 

O reproduces the up-up-down-down spin ordering. The method is free of the sign problem that is 

usually the main limiting factor for simulating fermionic systems in dimension higher than one. 
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I. INTRODUCTION 



Interacting electrons in crystals exhibit a variety of cooperative phenomena. This is 
particularly evident for Mott insulators whose valence electrons are localized due to a rather 
strong intra-orbital Coulomb repulsion. The degree of electronic localization is controlled 
by the ratio U /t, where t is the dominant transfer integral between different atomic orbit als. 
The strong {U/t ^ 1) and weak-localization {t/U <C 1) regimes can be treated by expanding 
in the small parameters t/U or U/t. Unfortunately, this is not true for the intermediate- 
coupling regime {U/t > 1) because of the lack of small parameter. Therefore, it is very 
important to develop controlled and unbiased numerical techniques that can address this 
elusive regime. The main obstacle for Monte Carlo simulations is the sign problem that 
can arise from the fermionic statistics of the electronic degrees of freedom or the frustrated 
nature of the Hamitonian terms. 

In this paper, we focus on the Vanadium spinels A^+V2^04~, a class of highly frustrated 
Mott insulators that belong to the intermediate-coupling regime. These systems comprise 
a pyrochlore lattice (corner sharing tetrahedra) of ions (see Fig. |l|(a)). Each ion 
has two in 3-d electrons in the t2g orbit als, while the and Zn^+ ions are in closed shell 
configurations. Consequently, the low-energy electronic spectrum consists of spin and orbital 
excitations that arise from the localized 3-d electrons in the V^^ ions. 

The lattice symmetry of AV2O4 is cubic (Frf3m) at high temperatures. The crystal field 
produced by the octahedral coordination of ions splits the 3rf-orbitals into high- 
energy eg and low-energy t2g orbitals. ZnV204 undergoes a structural cubic to tetragonal 
(Mi/amd) transition at T = 45 ^ 51(K) [K [2]. This structural transition causes further 
crystal field splitting of the t2g orbitals into x^-orbital and ^z, zx-orbitals. (See Fig{l](c).) 
A magnetic transition to up-up- down- down magnetic ordering (uudd-MO) occurs at a lower 
temperature T/v = 31 ^ 40(K) [1^ 2j. The magnetic ordering has a period of 4 lattice sites 
(• • • ttii • • •) ioT chains oriented along the yz (0, ±1, ±1) and zx (±1, 0, ±1) directions. The 
origin of the uudd-MO and the lack of orbital ordering in this material have been an open 
this magnetic ordering was originally reported in 1973 [3]. More recent measurements of the 
pressure dependence of T/v in diflFerent spinel vanadates [4J indicate that ZnV204 belongs to 
the intermediate-coupling regime: T/v decreases with pressure instead of increasing accord- 
ing to the Bloch's law [5j that is expected for the strong-coupling regime. Previous attempts 
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at explaining these properties were based on strong-coupling P-[9] or weak-coupling expan- 
sions [10] whose validity is not guaranteed for the intermediate-coupling regime relevant for 
ZnV204. 




FIG. 1: (a) up-up-down-down magnetic ordering at low temperature, (b) unit cell of pyrochlore 
lattice, (c) crystal field splitting due to structural transition. 



II. THREE-BAND HUBBARD MODEL ON THE PYROCHLORE LATTICE 



The three-band Hubbard model relevant for the family of spinel vandates AV2O4 is 
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where the operator rfj^^rr electron of spin a in the orbital a of the V^^ ion with coordinates 
'^5 ^olgT = dacrr^oiar^ ^ud CT^-^ = Yli ^ u d^niiV^livdaur ) whcrc cr"^ are the Pauli matrices, ta 
is the transfer integral between a-orbitals of nearest-neighbor ions, U is the intra- 
orbital Coulomb repulsion, J is Hund's coupling, ji^ is the chemical potential whose orbital 
dependence is a consequence of the crystal field splitting induced by the structural transition, 
l^yz = l^zx = l^xy — ^5 ^ud ^(j- r')^a Tcpresents summation over all pairs of the nearest- 
neighbor sites connected by a-type bonds. (As shown in Fig{l](b), the pyrochlore lattice has 



three types of bonds: xy, yz^ and zx.) This imphes that electrons occupying an a-orbital 
can only move along the single chain of a-type bonds that contains such orbital. Note that 
the on-site interaction term has only two independent parameters, U and J, because of the 
rotational symmetry of the single-ion Hamiltonian [llj . 

By expanding around the strong-coupling limit of Hsbh, Tsunetsugu and Motome derived 
an effective Kugel-Komuskii Hamiltonian whose ground state exhibits anti-ferro layered 
orbital ordering (00) |6j. Moreover, by adding an antiferromagnetic exchange between 
third-nearest-neighbor spins, they showed that an uudd-MO is stabilized and coexists with 
the 00 [7J. Since the 00 has not been observed in any member of the AV2O4 family, 
Tchernyshyov and Maitra-Valenti proposed alternative orbital orderings [8„, On the other 
hand, by starting from the opposite weak-coupling limit, Chern and Batista derived uudd- 
MO without any orbital ordering by using a mean-field theory p^. In this paper, we show 
that our unbiased approach not only explains the observed magnetic ordering but also the 
lack of the 00 in the intermediate coupling-regime relevant for the vanadium spinels. 

Quantum Monte-Carlo methods for the simulation of fermionic Hamiltonians such as 
^3BH typically suflFer from the well-known negative sign problem introduced by the fermionic 
statistics and by certain oflF-diagonal terms. To avoid the negative sign problem we replace 
Heisenberg-like Hund's interaction (fifth term in Eq. ([T])) by an Ising-like term because. In 
doing so, we are assuming that the optimal magnetic ordering of 'Hsbh is always collinear. 
In addition, we eliminate the hopping of singlet pairs between different orbitals of the same 
atom (sixth term in Eq. ([T])) by assuming that the probability of double occupancy of a given 
orbital remains rather low in the intermediate-coupling regime. The resulting Hamiltonian, 
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has only diagonal on-site interactions. Since all electrons are restricted to move along one- 
dimensional chains, there is no fermionic sign problem (fermions cannot be exchanged) 
and the only remaining challenge for efficient simulations is the the presence of geometric 
frustration. 
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III. METHOD 



We apply the world-line quantum Monte-Carlo method (QMC) [T2] to the Hamiltonian 
([2]). We use a modified directed-loop algorithm for the updating procedure [13]. In this algo- 
rithm, we insert a pair of discontinuities to the world-line configuration in (rf+l)-dimensional 
space, and move one of them stochastically by updating the world-line configuration. When 
the discontinuity comes back to the creation point where the other discontinuity is located, 
these discontinuities annihilate each other. Although the Hamiltonian ^ does not cause 
any negative sign problem, the geometrically frustrated nature of the underlying lattice gen- 
erates a "freezing problem" at low temperatures. To solve this problem, we introduce both 
thermal and quantum annealing in our QMC simulations. In the thermal annealing process, 
we start from an inverse temperature /3init and add 5/3 every A^ann Monte-Carlo sweeps. We 
choose /5init = 1.0(eV~'^), Aann = 300, and 5(3 = {(3 — Anit)/(Achange — 1) so that the inverse 
temperature is /3 after Achange temperature changes. Because of the rather strong interaction 
between electrons, the discontinuity is often geometrically trapped and localized. To avoid 
this "freezing of the simulation" , we apply a quantum annealing technique. We add a new 
term, 

+d^(xar^o.^rd^^cr'r^(^^'T ^ (3) 

to % to liberate the trapped discontinuity. We use the re-weighting technique to calculate 
the physical quantities of % from the simulation with Ti + Tiq. The procedure is rather 
simple. We sample the world-line configurations only without any vertices corresponding to 
T-Lq. We choose q SlS q/3 = 0.002 <C 1 so that we can sample 90% of world line configurations. 

We define the unit cell of the pyrochlore lattice as shown in Fig. [l](c). The primitive 
vectors are ai = (1/2,1/2,0), as = (1/2,0,1/2), and = (0,1/2,1/2). Each unit ceU is 
labeled hy R= n^a^ (n^ is integer). The coordinates of the four ions in the unit cell are 
ei = (1/4,0,0), 62 = (0,1/4,0), 63 = (0,0,1/4), and 64 = (1/4,1/4,1/4). To reveal the 
existence of uudd-MO and estimate the magnetic transition temperature, we calculate the 
susceptibility 
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^^R^ is the order parameter for uudd-MO, i.e., (^^) or (^^) finite in uudd-MO. 



IV. RESULTS 

We use [/ = 3.9 (eV), and J = 0.8 (eV) close to the previous theoretical works [6l l9l [T4]. 
[t is estimated as 0.35(eV) [I5].) Since we are focusing on the magnetic order transition, 
we assume that the crystal has tetragonal symmetry: A > 0. For concreteness, we use 
jjixy = 3.05 (eV) and A = 1.0 (eV) so that one electron occupies xy orbital and the other 
electron occupies the yz or zx orbitals at low temperatures. Figure [2] shows the results of 
our QMC simulations at t = 0.8 (eV). The finite size scaling of XsL~'^^^ has a crossing point 
at T = 0.145 ± 0.005 (eV) where the specific heat curves have a peak. Figure [3] shows five 
transition temperatures estimated by using the same finite-size scaling analysis. 




FIG. 2: Results of quantum Monte Carlo simulation with U = 3.9 (eV), J = 0.8 (eV), A = 1 
(eV), fixy = 3.05 (eV), t = 0.8 (eV), (a) specific heat , (b) finite-size scaling of Xs by assuming 
second order phase transition with 3D Z2 symmetry breaking (i.e., 77 =0.04). 
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FIG. 3: Phase diagram at = 3.9 (eV), J = 0.8 (eV) , A = 1 (eV), and /x^^/ = 3.05 (eV) with 
t — txy — tyz — tzx' Blue shaded regime is guide to hne. 

V. DISCUSSIONS AND CONCLUSIONS 

The transition temperature changes non-monotonically as a function of t. By calculating 
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where Arir = "^a^yzcrr — '^zxaV) we have also confirmed the 00 proposed by Tsunetsugu- 
Motome type does not appear for the range of parameters shown in Figure [3] This result is 
qualitatively consistent with the above mentioned experimental observations. On the other 
hand, the magnitude of T/v is much higher than the observed values, but it drops very rapidly 
for smaller values of t/U. The natural question is if the 00 is still absent for values of t/U 
such that T/v becomes comparable to the experimental values. An answer to this question 
requires of very low-temperature simulations for which the freezing problem becomes more 
challenging. It is important to note that the rather high value of A that we used for our 
simulations tends to increase T/v because of the localization of a single electron {S = 1/2) in 
the xy-orbital. We are also ignoring the vr-bond hopping of electrons that can reduce further 
the magnitude of T/v- Takubo et.al estimated the amplitude of tt hopping is half of t |15j . 

In summary, we obtained the observed uudd-MO without 00 for the intermediate- 
coupling regime of the three-band Hubbard model by using QMC simulations complemented 
by thermal and quantum annealing techniques. Preliminary results for t = 0.3 (eV) (not 



shown in this paper) are also confirming the existence of the 00 proposed by Tsunetsugu 
and Motome [6j. A complete the phase diagram that includes to the transition to the 
strong-coupling regime will be presented in a future work. 
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